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We reexamine the tunneling density of states in quantum dots and nanoparticles within the model 
which is extension of the universal Hamiltonian to the case of uniaxial anisotropic exchange. We 
derive the exact analytical result for the tunneling density of states in the case of arbitrary single¬ 
particle energy spectrum. We find that, similar to the case of the isotropic exchange, the tunneling 
density of states as a function of energy has the maximum due to a finite value of the total spin of 
the ground state near the Stoner instability. We demonstrate that there are no additional extrema 
which have been predicted on the basis of perturbative analysis [M.N. Kiselev and Y. Gefen, Phys. 
Rev. Lett. 96, 066805 (2006)]. 
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I. INTRODUCTION 

Quantum dots host rich physics that has been attract¬ 
ing experimental and theoretical interest for many years 
W- In the metallic regime when the Thouless energy 
(Axh) is much larger than the mean single-particle level 
spacing (5), E^\^/5 ^ 1, electrons in quantum dots can 
be universally described by an effective zero-dimensional 
Hamiltonian [^. The advantage of this so-called uni¬ 
versal Hamiltonian is the reduction of a set of matrix 
elements which describe electron-electron interaction in 
the single-particle basis to three parameters. They are 
the charging energy (Ac), the ferromagnetic Heisenberg 
exchange (J > 0) and the Cooper-channel interaction. 
Typically, the charging energy is large, E^, ^ S, and sup¬ 
presses a real electron tunneling through a quantum dot 
at low temperatures T ^ E^'. the phenomenon known as 
the Coulomb blockade [7]. It can be seen as suppression 
of the tunneling density of states at low energies |8l9j . 

Contrary to the charging energy, the exchange interac¬ 
tion typically is small, J ^ 5. Provided the size (L) of a 
quantum dot is larger than the Fermi wave length, the ex¬ 
change interaction can be estimated as J/S = —Fq where 
Fq denotes bulk value of the Fermi-liquid triplet channel 
interaction parameter. In bulk materials at Fq = — 1 
there is the Stoner instability corresponding to the quan¬ 
tum phase transition between a paramagnet and a fer- 
romagnet. Hence the ground state in quantum dots is 
expected to become fully spin polarized at J = 6. Sur¬ 
prisingly, an interesting regime with a finite total spin in 
the ground state is possible in quantum dots [S] . For the 
equidistant single-particle spectrum a partial spin polar¬ 
ization can be found at 6/2 < J < 5. With increase of 
J/S towards the Stoner instability, the total spin in the 
ground state increases monotonously. This phenomenon 
of subsequent transitions between the ground states with 
different values of the total spin is termed as the meso¬ 
scopic Stoner instability [3]. At J = 5 all electrons in 
a quantum dot become spin polarized. The mesoscopic 
Stoner instability is restricted to systems of a finite size 
and does not survive in the thermodynamic limit ^ > 0. 


The finite total spin of the ground state yields the 
Curie type behavior of the static spin susceptibility and 
its moments near the Stoner instability [TTVH5) . In¬ 
terestingly, existence of the mesoscopic Stoner instabil¬ 
ity can be seen in the low temperature electron trans¬ 
port through a quantum dot. The finite total spin in 
the ground state leads to an additional nonmonotonicity 
(in comparison with nonmonotonicities due to Coulomb 
blockade) of the energy dependence of the tunneling den¬ 
sity of states [TOUI^ and to enhancement of the shot noise 

[Sl¬ 
it is worthwhile to compare the case of isotropic 
(Heisenberg) exchange with the case of Ising exchange. 
For example, the latter can be realized in a two- 
dimensional quantum dot in the presence of strong spin- 
orbit coupling. A spin-orbit coupling invalidates the uni¬ 
versal Hamiltonian description of a quantum dot since 
fluctuations of the interaction matrix elements cannot 
be neglected even in the metallic regime, 6/ETh ‘C 1 
mm- In case of a two-dimensional quantum dot only 
in-plane components of the spin are mixed with the or¬ 
bital degrees of freedom whereas the perpendicular com¬ 
ponent of the total spin is conserved. If parameters of a 
quantum dot satisfy the following condition, {Xso/L)"^ ^ 
iETh/S){L/Xso)‘^ ^ 1, where Xso stands for a spin-orbit 
length, the low energy description can be given by the 
universal Hamiltonian with the Ising exchange {Jz > 0) 
mm- In this case there is no mesoscopic Stoner in¬ 
stability for the equidistant single-particle spectrum [3]. 
Since the total spin in the ground state is zero for all 
Jz < 6, the tunneling density of states is almost indepen¬ 
dent of Jz pn] . 

The simplest way to incorporate Heisenberg and Ising 
interactions is to consider the universal Hamiltonian with 
an uniaxial anisotropy of exchange interaction. Albeit 
this model is not fully microscopically justified it can be 
relevant for nanometer-scale ferromagnetic nanoparticles. 
We note that significant anisotropy of an exchange inter¬ 
action was revealed in experiments on tunneling spectra 
in such nanoparticles |2T]. The model which resembles 
the universal Hamiltonian with anisotropic exchange al- 
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lows to explain the main features of experimentally mea¬ 
sured excitation spectra |22j . Anisotropic exchange in¬ 
teraction in nanoparticles can be caused by bulk magne¬ 
tocrystalline, surface or shape anisotropy. The presence 
of spin-orbit scattering results in large mesoscopic fluctu¬ 
ations of the anisotropic part of the exchange interaction 
[23I24| . In quantum dots anisotropic exchange interac¬ 
tion can be induced by ferromagnetic leads pS] , 

The tunneling density of states for the universal Hamil¬ 
tonian with uniaxial anisotropy of exchange interaction 
was studied in Ref. m by means of the perturbation ex¬ 
pansion near the Ising case. It was found that anisotropic 
exchange interaction induces reentrant behavior (with 
two maxima and minimum) of the tunneling density of 
states as a function of energy. In contrast, energy depen¬ 
dence of the tunneling density of states is monotonous in 
the case of Ising exchange [20] and has the single maxi¬ 
mum in the case of Heisenberg exchange [mnj. Thus the 
result of Ref. uni indicates interesting physics due to fi¬ 
nite spin in the ground state in the presence of anisotropic 
exchange. However, this expectation does not supported 
by recent calculations of the spin susceptibilities for the 
model with anisotropic exchange m- The spin in the 
ground state monotonously reduces with change of the 
anisotropy from zero (in the case of Heisenberg exchange) 
to the maximal value (in the case of Ising exchange). 

In this paper we reexamine the problem of calculation 
of the tunneling density of states for the universal Hamil¬ 
tonian extended to the case of exchange interaction with 
uniaxial anisotropy. Within this model we derive exact 
analytical result for the tunneling density of states. We 
analyze this exact result in the cases of zero temperature 
and for temperatures larger than the mean single-particle 
level spacing. We demonstrate that 

(i) similarly to the case of isotropic exchange, in vicin¬ 
ity of the Stoner instability the tunneling density of 
states as a function of energy has the single maxi¬ 
mum due to presence of the finite total spin in the 
ground state; 


exchange interaction these corrections are negligible but 
results in rich physics beyond the zero-dimensional ap¬ 
proximation 1^. Also in the case of isotropic exchange 
interaction m the fluctuations of single-particle levels 
do not affect strongly the tunneling density of states. 

The outline of the paper is as follows. We start from 
definition of the model Hamiltonian and partial disen¬ 
tangling of spin and charge degrees of freedom (Sec. 0- 
In Sec. [ml we derive the exact analytical expression for 
the tunneling density of states. Zero temperature behav¬ 
ior of the tunneling density of states is explored in Sec. 
JVj The temperatures well above the mean single-particle 
level spacing are considered in Sec. jV] We conclude the 
paper with the summary (Sec. VI). 


II. FORMALISM 

A. Hamiltonian 

We consider the following model: 


H = Ho + Hc + Hs. (1) 

Here Hq describes noninteracting electrons, 

Hq = 'y ( y ( eadl^^Claa ■ ( 2 ) 

O' <T=± 


It involves the spin-independent single-particle energy 
levels Cq, and the single-particle creation (a)j^) and an¬ 
nihilation (oacr) operators. The next term He in the 
Hamiltonian Q takes into account the effect of direct 
Coulomb interaction among electrons in zero-dimensional 
approximation, 

Hc = Ec{n-Nof, h = ^ajj„.aa.o-, (3) 

Oi^G 


(ii) there are no additional extrema in the energy de¬ 
pendence of the tunneling density of states contrary 
to predictions of perturbative analysis of Ref. [T0| . 

In our analysis, we use the following standard simplifi¬ 
cations. We do not consider interaction in the Cooper 
channel which drives superconducting correlations in 
quantum dots [26] . This is allowed in the case of repulsive 
Cooper channel interaction [Bj. Although our exact ana¬ 
lytical result for the tunneling density of states in the case 
of uniaxial anisotropic exchange interaction is valid for 
an arbitrary single-particle spectrum, in its analysis we 
avoid consideration of randomness of single-particle lev¬ 
els. As we mentioned above in this case one needs to take 
into account corrections to the zero-dimensional Hamilto¬ 
nian which stem from fluctuations of the matrix elements 
of the electron-electron interaction in spite of the 

metallic regime, ^/i^Th ^1. In the case of isotropic 


where Ag denotes the background charge. The last term 
Hs describes ferromagnetic anisotropic exchange inter¬ 
action (J_L > 0, Jj > 0) 


Hs = -J^[Sl + SD - Jjl 



a,GG' 


Here cr = {cr^, cr*', cr^} stands for the standard Pauli ma¬ 
trices. The exchange part Hs of the Hamiltonian Q 
interpolates between the Heisenberg exchange, J_l = Jz 
and the Ising exchange, J_l = 0. In both cases, Eq. Q 
reduces to the universal Hamiltonian [6|. The Hamilto¬ 
nian 0 . with the Ising exchange is used for description 
of lateral quantum dots with strong spin-orbit coupling 

mm- 
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B. Partial disentangling of spin and charge 

Our aim is to compute the tunneling density of states 
^'(e) for the Hamiltonian 0. It can be conveniently ex¬ 
pressed via single-particle Green’s function in the Mat- 
subara time domain m 

oo 

i^(e) =cosh ^ j a,aa (it +13/2), (5) 

— OO ’ 

where /3 = 1/T. In the Lagrangian formalism, the Mat- 
subara Green’s function (matrix in the spin space) can 
be written as 

Ga(ri,r2) = J V[il, (j), 6 ]^a{Ti)^aiT 2 ) 

Z = [ (6) 


The Green’s function 5a(Ti2,0o) corresponds to the ac¬ 
tion Stot with (j) substituted by in the first line of 
Eq. Q and by 0 in the second line. Thus, Qa{Ti2,4’o) 
can be formally considered as the single-particle Green’s 
function for the Hamiltonian % = T-Lq + Hg where "Ho 
is given by Hq (see Eq. 0 ) in which is replaced by 
Eq.o- = Eq.o- ~ M + i4'o- We emphasize that the charge 
and spin degrees of freedom are not fully disentangled. 
The remnant trace of Hq is encoded in (j)Q which leads 
to a small imaginary shift of the chemical potential. We 
remind that the grand partition function for the Hamil¬ 
tonian Q is given as fTKm] 

^ ^ 2 ' ^-PE^(n-NQ)'^+Pnn+l3J^m(m+l) 

|m-|-l/2|-l/2 

X sgn(2m -I- 1) E (-^2) 

;=-|m-|-l/2|-|-l/2 


Here T denotes the time ordering and Stot is the imag¬ 
inary time action for the Hamiltonian 0 after the 
Hubbard-Stratonovich transformation: 


^tot = drlY, 


9r — 6 q ; + /)/ + + 


O’ • 0 




+ 


el + 01 01 


X I 

4J± 




(7) 


where = n/2 ± m, and is the Darwin-Fowler 
integral: 



III. EXACT EXPRESSION FOR THE 
TUNNELING DENSITY OF STATES 


We have introduced the chemical potential /i and the 
Grassmann variables = (V^at, '4’aiY', = (lAa^^ V’a;) 

to represent electrons on the dot. The scalar (j) and vector 
6 bosonic fields were used to decouple the direct Goulomb 
and exchange interactions, respectively. We start from 
performing a gauge transformation in the charging sector 
by splitting the field ^(t) as 


A. Wei-Norman-Kolokolov transformation 

In the Hamiltonian formalism Ga(Ti 2 ) can be written 
as 


1 j—ICa{—iT,—ir + iP), 

2 1 1Ca{-iT - i(3, -ir), 


r > 0, 
r < 0, 


(14) 


,, , 7 , , 27rTO 

0(r) = <j){T) + —— + (j)Q, 


dT(j){T) = 0, \(j)o\ < ttT (8) 


where Z = exp (—PH) and 


ICo 


,(4,t-) = Tre-*GH4 




(15) 


with integer m. The part ^(t) -|- 2'KmT of (/'(t) can be 
gauged away (see Refs. [ ^101201,^ for details). The 
Green’s function 0 becomes 

ttT 

Ga(Tl,r2)= J ^^^°E (ri2,4)ga(n2,^o), (9) 

—ttT 

ttT 

^ = (10) 

-ttT 


where ri 2 = ri — r 2 . The so-called Goulomb-boson prop¬ 
agator reads 


D{t, (j)o) 


^^-E^\T\{l-\T\T)+zcl,o{/3k+T)-l3E^(k-No+rT)^ ^ 
feGZ 

( 11 ) 


Using the commutativity of TLq and Hg we can split the 
evolution operator for "H into two parts, exp(IfH) = 
exp(ifHo) exp(itiJs). Next we apply the Hubbard- 
Stratonovich transformation to get rid of terms of the 
fourth order in electron operators in the exponent Hg: 


e 


= lim 

N^OO 


N 

[Y[der,]Y[Te^*^-^-/^ 


X exp 




n—1 


a 




(16) 


where A = t/N. Here and further on we omit nor¬ 
malization factors. They will be restored in the final 
result. In what follows we shall concentrate on evalu¬ 
ation of /CQ,(t+,t_). The corresponding partition func¬ 
tion Z has been already computed in Ref. [TS]. To 
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make further progress with Eq. ( |16[ ) we apply the Wei- 
Norman-Kolokolov transformation which allows 

us to rewrite T-exponent as a product of usual expo¬ 
nents: 


N 


Te 


iAOnSa _ 


exp is^A pp^r. 


N n 


exp tsP A ^ Kp P e 


n—1 


^-ipApp 


(17) 


n=l j = l 


where + ips^- Equation 0 is valid for both 

p = ±. We use the initial condition = 0. The 
variables 6 can be expressed via new variables Pp, and 


KpP as follows: 


^x,n '^P^y^n 


2 ~ ''‘p,ni '^z,n — Pp,n ''‘p,nV^p,n ' '>,n 

dx,n W0y,n _ Pp.n(^p,n + ^p,n-1) 

2 “ 2 




(kP -I- 12 

\^p,n ' ’^p.n—l) —r 


(18) 


A few comments are in order here. The vector in Eq. 
( [T7| ) is supposed to be real but the transformation ( [T8| 
suggests that it is complex. This corresponds to rotation 
of the contour of integration in Eq. 0- To preserve 
the number of independent variables we chose Pp^n to be 
purely imaginary, pp_„ = —p* and k+„ and to be 
complex conjugated, K+ = (kT"„)*. The transformation 


(181 assumes that the quantity {k^ ^ + >^p,N- _^)/2 corre- 


^p,N-l> 


spends to K^{t) in the continuous limit. In general, one 
can use any of discrete representations of K^(t) of the 
form vKppjq + (1 - ^ with 0 ^ z/ ^ 1. However, 

the symmetric one is special since for the choice z/ = 1/2 
it is sufficient to work with Eq. (16) to the first order in 


A. The Jacobian of the transformation (18) is given as 


exp(tpAX;^^iPp,n/2) [Ml- 

Rewriting two exponents in Eq. (15) with the help of 
representation 0 , we obtain 


ICo 


n n j dKl,n^dKp^^dpp^np exp 


p=± l,np = l 


ipApp^n^ 


X 1- 


Pp,i 


2J. 


^P.rip (^p.rip + ' 


''p.rip —1/ 


4J± 


'^P,np) \'^p,np + ^p,np-l) 




_ P’"P (,.P _ u-P \ 

^p,np —1/ 


n Tr[A(+)A(-)] 

7#a 


( 19 ) 


Here the limit Np —> c» is assumed. The quantity x = 1— 
J_l/ Jz characterizes a deviation from the case of isotropic 
exchange. 

The single-particle operators Aa ^ represent the evolu¬ 
tion operators. In accordance with Eqs. (0 - ^ they 
are defined as follows 


Np 

^(p) = exp I is^A^pp,p 


X exp 


Np 


is^A Y «^p,n exp(-tpA Y Pp,j] 

n=l j=l 


■ ( 20 ) 


Due to non-zero value of parameter x the action for k’s 


in Eq. (191 is not Gaussian: there are forth order terms. 


To get rid of such terms, we introduce a set of auxiliary 
variables Pp,np by employing the Hubbard-Stratonovich 
transformation: 


exp 


ipAx , \ 2 / „ p 

/^p,np} \^p,np + ^p,np-l) 


4J± 


ip Ax 
47r Jj_ 


X exp 


X J dpp^ripexp 

ip Ax 


ip Ax 2 


2Jj^ ^P,"p^p,rap (^p.rip 


^p,rip —1/ 


( 21 ) 


Next following Ref. [31], we introduce new variables 


K~P = 6“"- 

p,np Np,np^ 


<P = -yP 

P,np Np,np'' 


( 22 ) 


where 


Pp,np = -ipAxYiPp,n - Vp,n), 
n—1 

ip Ax 


(23) 


— A 


{Pp,np Vp,np)- 


Such choice of ap_„p and /3p,np allows us to remove terms 
of the second order in x’s and first order in p in the action 
in Eq. (19). It can be done within accuracy of the first 


order in A. We note that one needs to take into account 


the Jacobian of the transformation (23) 


Jp = exp 


-ipAx{pp^rip - »7p.np)/2 


(24) 


In the absence of magnetic field, = 0 and 

Aatt = Therefore, we concentrate on calculation 

of below. After evaluation of the single-particle 

traces in the expression (19) we find 
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N„ 


^an(4:i-)= III n / exp 

p—± tnp = l 


ipA ^ ^ 1 ipA[ 2 , 


-[(1 - xr)pp,„p + >f?7p,„p] - 


Pv,n„ 


4 r [f'P.np ^ _ 


X. 


-p 


- P'"p (yP - yP \ 

\'^p,rLp A.p,np —1/ 

A 


o-2iec,i + 




Afp 


L 2 


N„ 


N„ 


X cos 


ipA 

Z_^ PP:^P 


JVp 


PXp.ATp exp -ipA>f ^ (Pp,„^ - 77p,„J 


E +11® *''^+’’exp_ _ 

p—:tnp — l j p—:k Up — l 

N—p n—p 71—p 

+* A X-p.n_p exp -ipAx Y {p-p,n - V-p,n) + ip A Y P-P 

ri—p — l n=l n—1 

After integration over variables Xp,np (see details in Appendix B of Ref. [T2] ) we obtain 

P=± lnp = l J 7 


n„ = l 


(25) 


idz- 


J 2^1 = 1 


27rz2 


I g-u;„-2ze„t+ 


X ^ exp I —2^;^ cos 

p=± \ 


N„ 


?EE 


Pp,np 


p—it np —1 


Afp 


X I ^pe 

I p=± 


-ipAx E„p=i (Pp.pp -Pp.pp) 


AE 


p 'P’^pvx > e 

71„ — 1 


dpe ^exp|-fJ_L?;ay e 
-ipA J(l->f)pp,„ + ><Pp,p] 




(26) 


where 


= E 




—ie-y (t-i- —) 


75^a 


Wa = 




(27) 


75^a 


Now it is convenient to switch to continuous representa¬ 


tion. In order to transform expression (261 into a more 
standard form, we introduce new variabL 


es: 


ip{t) =ip [ dt'[{l - x)pp{t') + xpp{t')\ + Cp(0), (28) 

satisfying the following relations: 


E P “ ^ip(.tp) + ip^ / dtTjp{t) 


p=± 


= 0 , 


(29) 


E ^ = 0 - 

p=± 


Then after integration over variables rjp and introduction 
of auxiliary variable x we can write the functional integral 


for as the integral of the Feynman-Kac type: 


OO ^ 

/Can = / da:e-*+-+(‘+-‘+ [] | /I?[^p] 

E P=± o 


ygip lo^ d,tCp-{l-2ipx)^p{0)/2 I JJ" 

) l^OL 


\l2^| = l 


idz^ 

27rz2 


-y-Wc-2vp cosh[(^+(t+)-{_(t_))/2] 


J ^yva 

0 

x<5(C-H(i-H) +C-(/-) + 21n(4up)) 

X >cxj ^+ipxj^ / A)tp + [ip(tp)-(,p{0)]/2 

p=± 

Here the Lagrangians Cp are given as 

r - ^Lc2 _ +-L 

4 ■ 


(30) 


(31) 


Now it is more convenient to rewrite Eq. (30) in the 
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Hamiltonian representation: 

/ 

n 


_ ^ — 2ieat4 


7^a 



\|2:^ | —1 


p=± 


rfa; n I / d^pdC'p 


^—iJz>!x^ptp — {l — 2ipx)^'/2 


xS 1^5] Cp + 21n(4n„y)^ cosh[(c+-e 

x(e+|e-*«^‘+|e;)(e-|e*’'"‘ie-) 

Et 


-€-)/ 2 ] 


, if t — _^ ^PxJz^tp _j^xJ t 


P=± 


(32) 


corresponding to the Lagrangian (31) reads [34 


The Hamiltonian of one-dimensional quantum mechanics 

(33) 




nj = -J^—^+-e- . 

Its eigenvalues are given by Jiy^ and eigenfunctions are 
spanned by the modified Bessel functions K 2 ii/ where v 
is a real number: 

i^W) = sinh(27r^)it:2ji.(e“^'^^). (34) 

TT 

Using the following result (see formula 6.794.11 on p. 
743 of Ref. [55]) 

OO 

j dvvsmh{2^:v)K2^u{‘2va)K2^u{e~^+/'^)K2^u{e~^-/‘^) 

e+-e- 


0 

TT^ / 1 _i+±«^ 

= — exp-e 2 _ 2va cosh 

16 ^ V 4z;„ 2 


, (35) 


we can integrate over y, and Then we obtain 

[c = (c^-e;)/2] 


/ 




'l4a 


idz^ 1 e 


27rz2 I u, 


dxd(^ e 


C/2 


e t-) J (ii 2 K 2 i^{ 2 vo,) y ^ 

0 

X < g('^'^a — xxJz+ipxJz/i)tp^—ipJ±iyitp+ipJ±iy^t-p 


Xg2..pC-pC/2Q^^^(gPC/2) I 


(36) 


where 


OO 

Qu+,^.{z) = z f W_{i'p\^ + 2p\nz). (37) 

-OO P = ± 


Now we use the following identity (see formula 6.576.4 
on p. 676 of Ref. [35]) 

OO 

/* i^+A —1^1/ 

dxx~^K^{ax)K^{bx) = T {xpg) 


22+Ar(i_A) 11 

^ ^ p,q—± 

X 2Fi(r++,r_+, 1 - A; 1 - b'^/a^), 


(38) 


where Xpq = {1 — X + pp + qi')l2, r(a;) stands for the 
Gamma function, and 27^i(aj b, c; z) denotes the hyperge¬ 
ometric function. Then we obtain (t_|_ — t_ = —ij3) 




p —2zeat+ 




dxdC, e 


2ixQ—^Jz >cx^ 


dh sinh(/i) nna 

'yz^OL cr=± 


^—Be-y — ah 


p=± 


-(l+p)C/2 


^^^{ep,->^xJz+^pJz/A)tpy^(^2h + ipJxtp, C, ^ Jj_)] , (39) 

where the function W is defined as 

./X - 2ay\ 

i^JT) 


1 


4 sinh y 

sinhy 

4e“^ exp 

( (x- 2y) 
V iz 


erf ( 


■) 


(40) 


Here erf(z) = (2/-y7r) dt exp(—t^) stands for the error 
function. Also we restored all necessary normalization 
factors. Next integrating over x, C, h and (po in Eqs. 
© and ( [l4| ), we find the following result for the single¬ 
particle Green’s function: 


G.nir) = - E 


^—^Ec{n—Noy-\-l3pn+l3J±m{m-\-l) 


OO 

X.- ^ . p(J^-J^)r{l-rT) f dbe^^Al 

SZy/J, - J 

— OO 

y g-[eo-A«-l-Bc(2n-2Afo + l) + Jj_(m+l/4)-|-6/2]T _/g,S6/2 

xT{Pb,2m+ 1) Zn^{ea)Zn^^ - Zn^ + iZn^_i{ea) 

-T(-/36,-2m) Zn^Z„^{ea) - Zn^{ea)Zn^ |. (41) 

Here = n/2 ± m, and Zn(ea) is the integral of the 
Darwin-Fowler type: 


Zni^a) = J —e 

0 

and the function 


db p-i9n 
27r 


n (1 + e 




( 42 ) 


^ sinh(xz/2) 

sinh(z/2) x sinh^(z/2) ’ 
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Finally, substituting the expression (411 for the Green’s 
function into Eq. ([^ and performing integration, we ob¬ 


tain the following result for the tunneling density of states 
for the Hamiltonian ([^: 


/(e) = 


1 -I- e 




'Ziny'Zn 


„ —/3-Ec («—-^Vo) ^-I-/3«-I-/3 Jr m (m+1) 


|m-|-l/2|-l/2 


sgn(2m -I- 1) 


„/3(J.-Jr)F 


E 




i=-|m-|-l/2|-|-l/2 


m — I 
m 


X |^i5^e — ec(-|-/r — Ef.(2n — 27Vo -I- 1) — J±{'ni 1/4) -|- — J±)il + 1/4)^ 

-I- (5^e — Ca + ^ — Ec{2n — 2 Nq -I- 1) -I- JL(m -I- 3/4) -|- (J^ — Jl)(^ -I- 1/4)^ ^ 




Zn^ (2m 1)Z„ 

2m -|- 2 -|- 2 Z Zn^ (Ca) 1 


(44) 


This equation constitutes the main result of the present 
paper. We emphasize that our result is valid for an ar¬ 
bitrary single-particle spectrum {ca}. Each term in Eq. 
(441 corresponds to the tunneling of an electron with en¬ 


ergy e and a given spin into (from) a single-particle level 
with energy ea- Each delta-function describes the en¬ 
ergy conservation. The factor Zn(ea)/Zn measures the 
probability that the single particle level with energy Cq, 
is empty provided the total number of electrons is n. In 
the isotropic limit, Jz = J_l, Eq. (44) coincides with 
the result obtained in Refs. In the case of Ising 

exchange, J_l = 0, Eq. (44) transforms into the result 
of Ref. [10]. In the absence of exchange interaction, 
Jz = J± = 0 , the result (44) coincides with the expres¬ 
sion found in Ref. [9]. 

Due to breaking degeneracy of many-particle spectrum 
by exchange anisotropy, each delta-peak for isotropic case 
is replaced by 2m-I-1 peaks. Envelope of this set of peaks 
has the width of the order of 2m(Jz — J±)- As we will 
demonstrate below, this leads to smearing of the peak in 
the tunneling density of states in comparison with the 
case of isotropic exchange. 

By using the identities, ~ Zn(ea)] = nZn and 

Zn = Zn{ea) -f (cq,). One can check that Eq. 


(44) fulfills the sum rule: 


de 


^(g) 

1 -t 


= T 


dlnZ 


(45) 


IV. ZERO TEMPERATURE ANALYSIS 


We start analysis of Eq. (44) from the case of low 
temperatures T ^ S. For simplicity, we consider the 
case of Coulomb valley (Nq is close to an integer). Then, 
at r ^ 5 we can use the following relations: 


where E^^ stands for the ground state energy of n spin¬ 
less electrons and 0(a;) denotes the Heaviside step func¬ 
tion. In the case of the equidistant spectrum it is equal 
to En'^ = 5n(n — l)/2. Let us assume that the ground 
state of the Hamiltonian ([^ with Jz ^ J_l corresponds to 
the total spin S. Then in Eq. (44) one needs to take into 


account the contribution with m = S and I = ES only, 
ther contributions, e.g. with m = S—1 and I = E(S— 1), 
will be exponentially small. Hence, we find 

i/(e)= JzS - ^ 


1 

25--hi 


e„ - JzS - 


J± , Jz 


ea>e 


E -- 2 ■ 4 

5 

y] S[e^ + (2J^ - Jz)S 

E JzS+^+ 


J_L Jz 

+ — -h — 
2 4 


(47) 


ec,>e N„ 


Here ia = s + ^ — Ec — ea- It is convenient to rewrite 
as = ^ 5 ^ 1/2 - Es, where Es stands 

for the single-particle contribution to the energy of the 
ground spin with the spin S. Then let us introduce the 
following energies: 

£i=Es+i/2-Es-JzS-^-^, 

5 2 = As_i /2 -Es + JzS+'^-’^, 

5 3 = Es+1/2 -Es + JzS+^-^, 

£4 = — Es + (Jz — 2J±)S - y - -J. 

We note that the following inequalities: £3 > £ 1 , 2,4 and 
£1 ^ £ 4 , hold independent of the value of S. Therefore, 
only three different cases are possible: (a) £2 < £1 < £4 < 


(48) 




























FIG. 1: (Color online) Tunneling of an electron with spin up 
(left) and spin down (right) into a quantum dot with a finite 
value of spin in the ground state. 


v(i:)/vg 



FIG. 2: The sketch of dependence of the tunneling density of 
states on energy at zero temperature. The shaded areas are 
equal (see text). 


^3, (b) £i < £2 < £4 < £3, and (c) £i < £4 < £2 < £3- 
Which case is realized depends on the value of the total 
spin S in the ground state with = S. The energy 
£1 {£2) is the energy needed for an electron with spin 
up (down) which tunnels to the lowest available single¬ 
particle level (see Fig. [^. The energies £4 and £3 are 
required for the tunneling of a spin-down electron to the 
lowest single-particle level available for an electron with 
spin up. The energy £3 (£4) corresponds to the final 
excited state with the total spin S' — 1/2 (S -I- 1/2). 


As follows from Eq. (47) the tunneling is possible only 
if an electron energy e exceeds minjSi,£2,S3, £4}. Due 
to a finite value of the total spin in the ground state the 
tunneling of an electron is sensitive to its spin projec¬ 
tion. Initially, only electrons with one spin projection 
can tunnel into the dot. For electrons with very large en¬ 
ergies there is no dependence of tunneling on their spin 
projection. The characteristic energy that separates such 
large energies is given by £3. The sum rule (45) restricts 
possible behavior of the tunneling density of states. For 
the Coulomb valley and at low temperatures the sum rule 
(45) forces the integral J de i 2 {s) to be independent of the 
exchange interactions. As we shall see below this leads 
to the existence of the maximum in the tunneling density 
of states. In the cases (a) and (b) we find from Eq. (47): 


de 

£2 - £1 


£1 


^(g) 


0 , £1,2 < £ i 2 ) 

1 / 2 , £]^ ^ £1,2 < £12, 

< 1 , '^^12 ^ £1,2 < £4, ( 49 ) 

11 ^, £4 ^ £1,2 < £3, 

, 1 , £3 ^ £1.2, 


case (c) Eq. (47) yields 



^(e) 

1^0 


0, 

gl, 

2 < £l; 

! 


1/2, 

£1 


£1,2 

< 

£4, 

S-l-l 

2 S-I -1 ’ 

£4 


£^ 1,2 

< 


4 S-I -3 
4 S-I -2 ’ 

£2 


£^ 1,2 

< 

^ 3 , 

1, 

£3 


£^ 1 , 2 ' 




(50) 


Here vq = 2 S(s — ea) denotes the density of states in 
the absence of interactions. Energies £1^2 are measured 
with respect to Ec — /r. The sketch of dependence of the 
tunneling density of states on energy at T ^ 5 is shown 
in Fig. 1^ 

As follows from Eqs. (49) and (50) there is a step 
of height 1/2 between £1 and £2 in the envelope of the 
tunneling density of states for all three cases. How¬ 
ever, for the ground state with the total spin S « 
(J_L + Jz — 5)/[2{5 — Jz)] one can demonstrate that 
|£i — £2! < Jz, thus this 1/2 step will be smeared 
by temperature T > S. The maximum in the tun¬ 
neling density of states has the width of the order of 
£3 — £4 = J_l(2S' -I- 1) ^ ~ Jz)- Therefore, for 

temperatures in the range 6 T ^ j]_/{6 — Jz) one 
can expect that this maximum (of the relative height of 
1/(25') ~ [(5 — Jz]/J_L ^ 1) survives. We note that such 


temperature regime exists for ^/S{S — J^) J± only. 

We emphasize that the zero temperature analysis 
demonstrates clearly that the tunneling density of states 
has only single maximum. There is no other extrema in 
contrast to findings of Ref. uni based on perturbation 
theory in J±/Jz- 


V. THE TUNNELING DENSITY OF STATES AT 

t:$>5 


where £12 = min{£i,£2} and £12 = max{£i,£2}. For the 


Now we analyze the behavior of the tunneling den¬ 
sity of states a.t T 6. We neglect the effect of single- 
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particle level fluctuations. We restrict our analysis below 
to the case Jz ^ J±- We start from the integral rep¬ 
resentation (39) for the single-particle Green’s function. 
In order to perform integration over and h it is con¬ 
venient to express it through the noninteracting single¬ 
particle Green’s function. Then, provided the condition 
fj. ^ T ^ 6 holds we can approximate the noninteracting 


single-particle Green’s function as 

^ 1 _|_ g^(ea-At) S sinh(7rTr) 

Then, after integration over (j)Q, h we obtain 


(51) 


/(e) 1 -I- e 




1^0 


s,p—±: 


^ \ ^ ^ IQZs^yiJz - J±) 


db 


g-&V[4/3(^.-Ji)]eh/4 
cosh(&/4) sinh(6/2) 


^ sp¥ + , l(i{S-pJ±Y - j±) 


2’ 2l3{5-pJi_) ’V 4(5-jj_) ’ 


(52) 


where fl® = e + p — Ec{2n — 2 Nq + 1). Here we introduce 
the following functions: 


E{x,y) = 

and 

¥{x,y,z,u) = 

The quantity 


dt 




2 cosh(7rt) 


dt e’^(*‘+i/2) 
2 cosh(7rt) 


(53) 


erfi(z(2/-*t))e“'(‘'+i/4). 

(54) 


Zc = y^ 


,-PE^{n-No)‘^+PlJ.n 


(55) 


denotes the contribution to the grand canonical parti¬ 
tion function in the absence of exchange interaction and 
single-particle levels. The contribution due to exchange 
interaction is given by m 


Zs = 


S-Jz 


1/2 


ISJ1 


g4(6-j^) jr I 


{S-JxV 


//3J*), (56) 


where J, = (5 — J±){Jz — J±)/{S — Jz) and 

OO 


(57) 


Using the following asymptotes of the function J-{x, y): 


J"ix,y) = X + o 


1 1 j tanh y ^ 1, 


2 2 1 erf(l^), y » 1, 


(58) 


and the following relations J-{—x,y) = e ^T{x,y) and 

dy¥{x,y, z,u) = (y+i/2) jz(^x — 2yz^, \Jz"^ — u"^) 

Vtt 

( 59 ) 


one can check that in the case of isotropic exchange, 
Jz = Jx, the result ( [^ coincides with the result ob¬ 
tained in Ref. [T^]. As follows from Eqs. (58) and 
(59) the functions T{x,y) and ¥{x,y, z,u) deviates only 
slightly from the Fermi function fpix) = 1/(1 + exp(a;)) 
for y, z,u 1. Although u is always small, the parame¬ 
ters y and z can be large for J± close to S at temperatures 
(5 +J_l)V[4(<5-^±)>U»5. 


The plot of the energy dependence of the tunnel¬ 
ing density of states is shown in Fig. In accor¬ 

dance with zero temperature analysis there is the max¬ 
imum of ;/(e) at non-zero values of Jl and temperature 
T {6 + J_l) 2 /[ 4(5 — Jj_). We note that for the solid 
curve in Fig. we choose precisely the same parame¬ 
ters as in Fig. 2 of Ref. [101 ■ There are no additional 
extrema contrary to conclusions of Ref. uni based on 
perturbation theory in J±/Jz- 


VI. CONCLUSIONS 


In this paper we studied the Hamiltonian which is an 
extension of the universal Hamiltonian to the case of 
uniaxial anisotropic exchange interaction. Within this 
model we have derived exact analytic expression (see Eq. 
(El)) for the tunneling density of states for arbitrary 
single-particle spectrum. For {6 — Jz) ^ J± < Jz we 
analyzed the energy dependence of the tunneling density 
of states for the equidistant single-particle levels at low 
(T <C 6) and high {T 3> S) temperatures. In both cases 
we demonstrated that in addition to non-monotonicities 
due to Coulomb blockade there is the maximum in the 
tunneling density of states at the characteristic energy of 
the order of J'j_/{5 — Jz) (see Fig. |^. The relative height 
of the maximum scales as (5 — J^j Jx- Qualitatively, in 
the case of anisotropic exchange the tunneling density of 
states has the same energy dependence as for the case of 
isotropic exchange. Our findings demonstrate that ad- 
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FIG. 3: The tunneling density of states in the Coulomb valley. 
The solid (dashed) curve corresponds to Jz/5 = 0.92, J±/5 = 
0.85 and TjS — 1.02 (T/d = 3). The dotted curve corresponds 
to Jz/5 = J^/5 = 0.92 and T/5 = 1.02. 

ditional extrema in energy dependence of the tunneling 
density of states related with the anisotropic exchange 
interaction obtained in Ref. HO] on the basis of pertur¬ 
bative analysis in J^_/ Jz do not exist. 


The most promising regime for experimental investi¬ 
gation of non-monotonicity due to exchange interaction 
in the tunneling density of states is vicinity of the Stoner 
instability. For example, one can perform the scanning 
tunneling microscopy in system of nanoparticles made of 
nearly ferromagnetic materials, e.g. Pd or Pt with Co or 
Ni impurities, various transition-metal alloys with dis¬ 
solved Fe or Mn atoms, and rare-earth materials [3S]. A 
very promising candidate could be nanoparticles made 
from YFe 2 Zn 2 o compound which has the exchange inter¬ 
action J « 0.94(5 m- For quantum dots fabricated in 
two-dimensional electron systems exchange interaction is 
typically not large, J < S/2 and, therefore, can be impor¬ 
tant only at T < 5 [38]. However, recent experiments on 
two-dimensional strongly interacting electron systems in 
Si-MOSFET revealed the existence of electron droplets 
with finite spin of the order of 2 at low temperatures 
and low densities [39]. The well-known enhancement of 
electron-electron interaction in the triplet channel [10] 
suggests that the physics of mesoscopic Stoner instabil¬ 
ity in disordered electron system [41] can be relevant for 
Si-MOSFET at low electron densities [ 22 ] • 
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